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Abstract 

We report on the first far-IR detection of H^D" 1 ", using the Infrared Space Observatory, in the line 
of sight toward Sgr B2 in the galactic center. The transition at A=126.853 (im connecting the ground 
level of 0-H2D+, li.i, with the the 2%^ level at 113 K, is observed in absorption against the continuum 
emission of the cold dust of the source. The line is broad, with a total absorption covering 350 km s _1 , 
i.e., similar to that observed in the fundamental transitions of H2O, OH and CH at ~179, 119 and 
149 /im respectively. For the physical conditions of the different absorbing clouds the H2D+ column 
density ranges from 2 to 5xl0 13 cm -2 , i.e., near an order of magnitude below the upper limits obtained 
from ground based submillimeter telescopes. The derived H2D + abundance is of a few 10~ 10 , which 
agrees with chemical models predictions for a gas at a kinetic temperature of ~20 K. 
Subject headings: Astrochemistry — molecular processes — line : identification — ISM : molecules 
- ISM : individual (Sgr B2) — infrared : ISM 



1. INTRODUCTION 

The H^" molecular ion is a key molecule for the gas 
phase chemistry in interstellar clouds (Herbst & Kem- 
plerer 1973; Watson 1973). Produced by the fast reaction 
of H2 with H^, it reacts with almost all neutral atoms 
and molecules, and thus it is the precursor of a large 
number of complex neutral and ionic molecular species. 
Initially suggested to be present in molecular clouds by 
Martin et al. (1961), there is a long history behind the 
search for this crucial molecule. Due to its lack of perma- 
nent dipole moment, H^" must be observed in the near-IR 
through its ro- vibrational transitions (Oka 1980). The 
first searches only provided upper limits to the col- 
umn density (Geballe & Oka 1989; Black et al., 1990). 
However, H^" has been finally detected toward GL2136 
and the W33A clouds (Geballe k, Oka 1996). In addition, 
H^" has also been detected toward the diffuse interstellar 
medium (ISM) by McCall et al. (1998), and more re- 
cently in the line of sight toward several galactic center 
clouds (Oka et al. 2005). Unfortunately, these obser- 
vations only provide information on the H;j~ abundance 
in the most diffuse and warm external layers of molecu- 
lar clouds, and are limited to lines of sight with bright 
near-IR background sources. 

A different approach to study H^ in denser regions 
obscured by dust extinction is to observe the asymmet- 
rical species H2D + and HD^ , which do have a rotational 

1 Based on observations with ISO, an ESA project with instru- 
ments funded by ESA Member States (especially the PI countries: 
France, Germany, the Netherlands and the United Kingdom) and 
with participation of ISAS and NASA. 

2 Department of Physics, University of Lethbridge, 4401 Univer- 
sity Drive, Lethbridge, Alberta, T1J 1B1, Canada 



spectrum. These species play a crucial role in the deuter- 
ation enhancement of many molecular species in clouds 
with kinetic temperatures below ~30 K (see, e.g., Ger- 
lich et al., 2002 and references therein). In particular, 
the H2D + /H^" abundance ratio is expected to dramati- 
cally increase above the elemental D/H ratio as the gas 
cools down. The first calculation of the H2D 4 " rotational 
spectrum was reported by Dalgarno et al. (1973), who 
suggested that it could be detectable in dense and cold 
clouds. The H2D+ I1.0-I1.1 transition has been observed 
in the laboratory at 372.42134 GHz (Bogey et al. 1984; 
Warner et al. 1984). This precise measurement opened 
the quest for the detection of H2D + in dense molecular 
clouds. Initial searches for the o-H2D + 11,0-11,1 line, and 
also the p-H 2 D+ lo,i-0o,o (at 1370 GHz), performed with 
the Kuiper Airborne Observatory (Phillips et al., 1985; 
Pagani et al., 1992a; Boreiko & Betz, 1993), and with the 
Caltech Submillimeter Observatory (van Dishoeck et al., 
1992) were, however, unsuccessful. Fortunately, H2D+ 
has finally been detected toward the young low-mass pro- 
tostar NGC1333-IRAS4A (Stark et al. 1999), and in sev- 
eral dark clouds (Caselli et al. 2003; Vastel et al., 2004, 
2006a,b; Hogerheijde et al., 2006; Harju et al., 2006), 
and tentatively in protoplanetary disks (Ceccarelli et al., 
2004; see also Guilloteau et al., 2006). 

Observational evidence for the crucial role of H2D+, 
HD^ (Vastel et al., 2004), and even of D^, in the deuter- 
ation processes of the ISM have been widely probed by 
the detection of double and triple deuterated species such 
as D 2 CO (Turner, 1990; Ceccarelli et al., 1998), NHD 2 
(Roueff et al., 2000; Loinard et al., 2001; Gerin et al. 
2006), ND 3 (van der Tak et al., 2002; Lis et al., 2002), 
CHD 2 OH (Parise et al., 2002), CD 3 OH (Parise et al., 
2004), D 2 CS (Marcelino et al., 2005). 
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Ground-based observations of H2D + and HDj have al- 
ways been performed toward very dense and cold cores 
as the high Einstein coefficients of the involved transi- 
tions require high volume densities to produce significant 
emission. Therefore, it is difficult to estimate the H 2 D+ 
abundance in lower density regions through submillimc- 
ter observations. In this Letter, we report on the first 
detection of H 2 D + in the far-IR spectrum of Sgr B2 ob- 
served with the Infrared Space Observatory (ISO). The 
observed transition connects the ground state, li ; i, with 
the 2i 2 level at 113.4 K, and has a wavelength of 126.853 
/im. This detection opens the possibility to detect H2D + 
in clouds with moderate volume density in which the lo,i- 
li,i transition at 372 GHz will be subthermally excited 
and its intensity too weak to be detected with currently 
available radio telescopes. It also opens the quest to 
detect H2D+ toward far-IR luminous galaxies and high 
rcdshift objects. 

2. OBSERVATIONS 

The data presented in this paper are based on the fi- 
nal calibration analysis of the high spectral resolution 
line survey of Sgr B2 obtained with the Long Wave- 
length Fabry-Perot Spectrometer (LWS/FP) on board 
ISO in its L03 mode (Polchampton 2002; Polchamp- 
ton et al., 2007). The data for NH2, which are neces- 
sary to interpret the H^D" 1 " absorption reported in this 
paper, were previously presented by Goicoechea et al. 
(2004). The spectral resolution of 'the LWS/FP spec- 
trometer is ~ 35 km s _1 . The telescope was pointed to- 
ward a(2000)=17 /l 47 m 21.75 s ; <S(2000)=-28° 23' 14.1". 
The H2D" 1 " and NH2 data are presented in Figs. 1 and 2. 

The procedure to reduce the LWS/FP data have been 
described in detail by Polchampton et al. (2007). This 
provides a significant improvement over the previous re- 
duction of Goicoechea et al. (2004), who did not apply 
any interactive processing to the standard pipeline data 
products except for the removal of glitches, averaging 
of individual scans and removal of baseline polynomials. 
The improvements applied in the full spectral survey re- 
duction include determination of accurate detector dark 
currents (including stray light), improved instrumental 
response calibration (including a correction to account 
for adjacent FP orders) and interactive shifting of indi- 
vidual miniscans (to correct for uncertainty in the angle 
of the LWS grating). The signal-to-noise ratio was also 
enhanced by including every available observation in the 
survey (in this case 2 independent observations). The 
broad absorption feature shown in Figure 1 was previ- 
ously assigned only to p-NH2 (Goicoechea et al., 2004; 
transition 22,1-11,0 J=5/2-3/2). However, the final re- 
duction of the LWS/FP survey toward Sgr B2 shows that 
the broadening is real. The feature appears in 2 inde- 
pendent observations : ISO TDT number 50601112, ob- 
served on 1997 April 5 with LWS detector LW2; and ISO 
TDT number 50800515, observed on 1997 April 7 with 
LWS detector LW3. The fact that the feature appears on 
2 separate detectors rules out any spurious instrumental 
broadening of the p-NH2 line. 

Taking into account the latest data reduction, the 
broad feature detected at ^126 /im cannot be due to 
P-NH2 alone. In fact, NH2 far-IR lines are narrow and 
appear only at the velocity of Sgr B2 (see, van Dishoeck 
et al., 1993, for the detection of NH 2 in the submm). 
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Fig. 1. — Observed ISO spectrum at 126.853 (im toward Sgr B2. 
The upper panel shows the data and the expected contribution of 
P-NH2 (see Figure 2). The middle panel shows the 0-H2D+ line 
profile after subtracting the contribution from P-NH2. Finally, the 
lower panel shows the absorption of CH for comparison purposes. 

Only a few lines remain unidentified in the ISO far-IR 
line survey of SgrB2. The feature discussed in this pa- 
per must correspond to a line arising from a low energy 
level or to a blend of lines. However, molecules heavier 
than H 2 D + could have other lines in the far-IR domain. 
We have considered the possibility of b-type transitions 
from slightly asymmetric rotors like HNCO and HOCO + . 
Both have one transition at similar wavelengths. How- 
ever, the line of HNCO involves two high energy levels 
and can be easily ruled out. HOCO + has several other 
lines in the 120-127 /im range involving low energy lev- 
els. None of which are detected in the far-IR line survey 
of Sgr B2 (Polchampton ct al., 2007). Hence, we are 
confident that the broad feature detected at 126.853 /im 
corresponds to H 2 D + . 

The intensity scale in figures 1 and 2 was obtained 
by dividing the observed flux by a fitted continuum as 
described by Polchampton ct al. (2007). Hence, it is 
directly related with the line opacity convolved with the 
instrumental spectral profile of the LWS/FP spectrome- 
ter. The velocity scale in Figure 1 has been obtained as- 
suming A res4 = 126.853 /im (2363.325 GHz), i.e., the wave- 
length of the o-H 2 D+ 2i 2 -lii transition (Amano & Hirao, 
2005). The expected error on the frequency of this tran- 
sition is ~ 5 MHz, i.e., a velocity uncertainty < 2 km s _1 
(3cr) , which is much lower than the spectral resolution of 
the LWS/FP spectrometer. The p-NH 2 2 2 ,i-li, J=5/2- 
3/2 line is separated by ~-123 km s _1 with respect to 
the 2i2-ln line of o-H 2 D+. 
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3. RESULTS AND DISCUSSION 

The upper panel of Figure 2 shows the three spin- 
rotational components of the 0-NH2 22.0-I1.1 transition 
at ~117 /Ltm (Goicoechea et al. 2004). In order to com- 
pute the expected contribution of the P-NH2 22,1-11,0 
triplet at ~126 /jm we have modeled the NH2 absorp- 
tion assuming a rotational temperature of 20 K, an in- 
trinsic linewidth of 10 km s _1 , a total NH2 column den- 
sity of 1.2 xlO 15 cm~ 2 (Goicoechea et al. 2004), and 
an ortho/para ratio of 3. NH2 wavelengths and line 
strengths are from the Cologne Database for Molecu- 
lar Spectroscopy (Muller et al. 2001). The NH 2 hyper- 
fine structure is unresolved by the LWS/FP. The syn- 
thetic spectrum, convolved with the LWS/FP spectral 
resolution, is shown in Figure 2. The expected contribu- 
tion from P-NH2 to the absorption feature is shown as 
a continuous line in the upper panel of Figure 1. The 
middle panel of this Figure shows the same data with 
the expected contribution from P-NH2 removed. The re- 
maining absorption feature is broad and has several ab- 
sorption velocity peaks, which are similar to those found 
in unsaturated ground-state transitions of light species 
such as the (2,3/2)-(l,l/2) line of CH (see bottom panel 
of Figure 1). Therefore, in the following, we assign this 
unidentified broad feature to the o-H2D + 2i 2-li,i line. 

H2D+ is an asymmetric molecule with /j> a ~0.6 D (Dal- 
garno et al., 1973). It has two different species, ortho 
(Ka=odd), and para (Ka=even), with a spin weight ra- 
tio of 3:1. The three lowest energy levels, li,i, li.o and 
2i.2, which are at 0, 17.9 and 113.4 K (see insert in the 
middle panel of Figure 1), are enough to compute the 
partition function (F par ) at low temperatures. Hence, 
F par (T r ) ~ 3(l+e~ 17 - 8 / Tr ) is a good approximation for 
T r < 30 K, where T r is the rotational temperature. For 
example, F par (3 K)= 3.004 (all molecules in the ground 
state), while F par (30 K)=4.8. Although the li, level 
starts to be populated in this case, the 2i, 2 one is still un- 
populated. Consequently, if there is a strong continuum 
source behind or inside the source, we could expect to 
see the 2i,2-li,i transition of 0-H2D+ in absorption. The 
same behavior is observed in CH toward Sgr B2, where 
only the ground-state line is observed with the line pro- 
file shown in the bottom panel of Figure 1 (Goicoechea 
et al. 2004; Polehampton et al. 2007). 

The Einstein coefficient of the o-H2D + 2i.2-li,i tran- 
sition is large, 1.7xl0~ 2 s — 1 , which implies that the ro- 
tational temperature of the line will be extremely low 
in most relevant cases (or close to the dust temper- 
ature if continuum emission is optically thick), even 
when the I1.0-I1.1 transition is already collisionally ex- 
cited (its Einstein coefficient is two orders of magni- 
tude lower). Therefore, in order to derive a column 
density from absorption measurements of the 21,2-11,1 
line, we can assume that all molecules are in the li,i 
ground level, and that the transition has a low ro- 
tational temperature in absence of infrared pumping. 
The opacity of the 2i.2-li,i line can be written as 
r(T r ) ~ 0.067(N/10 13 )(10/Av)/F par (T r ), where Av is 
the linewidth at half intensity in km s^ 1 , and N is the 
column density of o-H 2 D+ in cm" 2 . If N(o-H 2 D+) = 10 13 
cm" 2 and Av=10 km s -1 , r(2 1)2 -l 1) i) would be 0.024 for 
T r =3 K and 0.014 for T r =3G K respectively. For com- 
parison purposes, the opacity of the li o-li 1 transition 
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Fig. 2.— Rest ISO/LWS/FP spectrum of Sgr B2 around 117 and 
126 (an, corrected by the averaged velocity of far-IR lines toward 
the source (62.7 km s _1 ). Lines due to 0-NH2 (upper panel) and 
P-NH2 plus o-H2D + (lower panel) are shown. The expected NH2 
absorption is shown by continuous lines in both panels (see text). 
An NH2 ortho/para ratio of 3 has been assumed. Note that the 
NH2 lines show absorption only at the velocity of Sgr B2. 

at 372 GHz would be 0.024 and 0.006 for T r =3 and 30 K 
respectively. 

The interpretation of the absorption produced by o- 
H2D+ in Sgr B2 requires knowledge of the rotational 
temperature of the transition. The cloud has a large dust 
opacity in the far-IR, Td(100/mi)=3.5 (Cernicharo et al. 
1997, 2006; Goicoechea et al. 2004). Hence, infrared 
photons cloud play an important role in the pumping of 
the molecular levels of 0-H2D+. In addition, the fraction 
of the cloud depth that is " seen" at this wavelength cor- 
responds only to its external layers, where th(127 /xm)~l, 
i.e., N(H 2 )~3xl0 23 cm" 2 , rather than the >10 24 cm" 2 
corresponding to the total column density along the line 
of sight (7^(127 /Ltm)~3). In order to estimate the ra- 
diative and collisional effects on the excitation of the 
2i,2-li,i transition, we have assumed de-excitation col- 
lisional rates £7(li,o-li,i)=c(2i,2-li,o)= 10~ 10 cm 3 s" 1 , 
and cr(2i,2-li,i)=5xl0~ n cm 3 s _1 . For a gas at 30 
K purely excited by collisions with N(o-H2D + )= 10 13 
cm~ 2 , and Av=10 km s _1 , the rotational temperature 
varies from 9 to 12 K for H2 densities of 10 5 and 10 6 
cm -3 reS pectively. Reducing the collisional rates by a 
factor 10 produces rotational temperatures of 6 and 8 K 
for the same volume densities. 

The role of radiative pumping by far-IR dust pho- 
tons can be estimated assuming that the gas we observe 
in absorption, Td(127/xm) = l, surrounds a region with 
T d (127/mi)=2.0 and T d =30 K. If the radius of the ab- 
sorbing shell is twice that of the region emitting the bulk 
of the infrared emission, R a h s /R cont =2, we obtain an ex- 
citation temperature of 17 K, with a little dependency 
on the assumed collisional rates. For R a f, s /R co „ t =3 and 
1.5, the rotational temperature is 15 and 20 K respec- 
tively. In the following we assume T r =10-20 K, and 
we use these values to derive upper and lower limits 
to the opacity of the H2D 4 " 2i,2-li,i line. Assuming a 
linewidth of 10 km s _1 we obtain r(2i 2-I1 1)=0.18 and 
0.36 for T r =10 and 20 K respectively. The correspond- 
ing column densities are ~ 9xl0 13 and 2.3 xlO 14 cm~ 2 , 
and the abundance of o-H2D + , for a total gas column 
density of 3xl0 23 cm" 2 , is ~ 3xl0" 10 for T r =10 K 
and 9xl0~ 10 for T r =20 K. The latter can be ruled out 
from the upper limit obtained by Pagani et al (1992a) 
for the emission of the I1.0-I1.1 line at 372 GHz, N(o- 
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H2D + )<l.lxl0 14 cm -2 . Hence, a low rotational tem- 
perature, T r ~10 K, and an abundance of 3xl0~ 10 are 
the best fit to the o-H 2 D+ absorption at the velocity of 
Sgr B2. 

For the other absorbing clouds which are not directly 
associated to Sgr B2, it is reasonable to assume radiative 
thermalization with the cosmic background, T r ~3 K, 
i.e., no significant collisional excitation or infrared pump- 
ing. Assuming an intrinsic line width of 10 km s" 1 , the 
observed absorption features at -100, 10 and 120 km s" 1 
have opacities of of 0.1, 0.15 and 0.06 and column densi- 
ties of ~ 4.5, 6.7, and 2.7 10 13 cm" 2 respectively. These 
clouds have low volume densities and total gas column 
densities of a few 10 22 -10 23 cm" 2 . The associated o- 
H 2 D + abundances are of a few 10 _10 -10~ 9 . 

The derived abundance of o-H 2 D + in Sgr B2 corre- 
sponds to the predictions of chemical models for Tg 
<20 K (see Pagani et al., 1992b; Walmsley et al., 2004; 
Flower et al., 2006). These models predict that the para 
species could be more abundant than the ortho one for 
that temperature and the expected density of the gas (a 
few 10 4 -10 5 cm" 3 ). Hence, the fundamental line of p- 
H2D + could produce significant absorption at 218.8 /xm 
in the direction of Sgr B2. The abundances observed 
toward the other clouds in the line of sight of Sgr B2 
indicate lower kinetic temperatures, Tk —10 K. In these 
clouds the predicted ortho/para abundance ratio favors 
the ortho species. Due the high dependency of the or- 
tho/para ratio with the gas temperature, simultaneous 
observation of both the 126.853 (im (o-H 2 D+) and the 
218.8 /im (p-H 2 D+) lines will provide a direct measure- 
ment of the efficiency of the reactions of H 2 D + with H 2 
and of the conversion of ortho-to-para H 2 D + and vice- 
versa, i.e., to check chemical models. 

Due to the poor spectral and angular resolution of ISO, 
and due to the limited sensitivity of the LWS/FP in- 
strument, the o-H 2 D + line at 126.853 /mi could only be 
expected in sources with a large column density of cold 
gas with broad linewidths in the line of sight. As an 
example, a broad absorption feature has been also re- 
ported in the bright IR galaxy Arp220 at 127 fim by 
Gonzalez- Alfonso et al. (2004) using the LWS grating 
spectrometer (spectral resolution more than 20 times 
lower than the LWS/FP). These authors assigned it to 
NH 3 . While several lines of NH and NH 3 have been ob- 



served in this source, no lines have been assigned to NH 2 
so far. Hence, the broad absorption feature assigned to 
NH3 by Gonzalez- Alfonso et al. , could also contain some 
contribution from o-H 2 D + . The detection of H 2 D + pre- 
sented in this work, shows that far-IR ortho- and para- 
H 2 D + lines can be used to trace the deuteration fraction- 
ation of Hjj~ in galaxies (see also Ceccarelli and Dominik, 
2006). Therefore, far-IR observations with the PACS 
instrument on board Hcrschel satellite (or in the future 
with SOFIA), will be an important tool to understand 
the cold ISM chemistry in distant galaxies where the D/H 
ratio will be different and the gas temperature affected 
by the local cosmic background temperature. 

Finally, the opacity of the 2i i 2-li,i transition for a dark 
cloud with Av=0.5 km s" 1 (see Vastel et al., 2004), N(o- 
H 2 D+) = 10 12 cm" 2 , and T r = 3 K is -0.05. While the 
expected emission in the 11,0-11,1 transition at 372 GHz 
will be negligible, the expected absorption at 126.853 
/j,m will be 5%, i.e., easy to detect with the expected 
sensitivities of the high spectral resolution instruments 
planned for SOFIA. Dark clouds (see, e.g., Nisini et 
al., 1999), prestcllar cores, and protoplanetary disks (see 
Wyatt et al., 2003) have enough continuum emission in 
the far-IR, hence, the 2i, 2 -li,i transition can be used 
MS ("I key observing tool to trace column densities of o- 
H 2 D + larger than a few 10 11 cm -2 , and to understand 
the deuterium fractionation at large spatial scales in cold 
dark clouds. Observations of the fundamental lines of p- 
HD+ at 691.6605 GHz and 0-HD+ 203.03 fim will provide 
additional constraints on these deuteration processes in 
molecular clouds and protoplanetary disks. 
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